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Introduction 


Accurate  prediction  of  sound  absorption  in  sea  water  is  essential  to  the 
effective  design  and  operation  of  sonar  systems  end  has  been  recognized 
as  a  serious  problem  for  many  years.  By  the  end  of  WWII,  it  hod  become 
clear  that  absorption  at  sonar  frequencies  is  an  order  of  magnitude  greater 
than  that  in  fresh  water.  Subsequent  propagation  experiments  provided 
accurate  estimates  of  the  magnitude.  By  the  early  1950's,  laboratory 
resonator  experiments  had  identified  the  cause  as  a  magnesium  sulfate 
relaxation  and  details  of  the  mechanism  and  relaxation  parameters  hod 
been  worked  out  [1,21.  In  1962,  Schulkin  and  Marsh  proposed  a  formula  for 
sea  water  absorption  based  on  field  experiments  in  the  frequency  range 
2-30  kHz  and  the  measured  relaxation  parameters  [3j. 

In  1965,  Thorp  {4]  reported  sound-channel  propagation  experiments  in 
the  Bermuda-Eleuthera  area  indicating  an  anomaly  at  lower  frequencies. 
The  excess  loss  was  fitted  by  addition  of  a  1  kHz  component  to  the  S&M 
formula.  The  result  became  known  os  the  "Thorp  formula"  (5). 

Mediterranean  experiments  were  reported  by  Leroy  (61,  showing  a  similar 
anomaly  but  with  somewhat  higher  magnitude  and  relaxation  frequency. 
Experiments  carried  over  the  next  two  decades  in  other  areas  confirmed 
the  high  degree  of  variability  of  the  extra  loss  and  regional  dependence  has 
therefore  become  much  more  critical.  The  major  factor  involved  in  this 
variability  has  been  identified  as  pH  (7). 

The  principal  chemical  relaxation  responsible  for  the  pH-dependent  loss 
has  been  shown  to  involve  boric  acid  (8).  The  chemical  mechanism  has  been 
identified  os  the  boric  ecid/carbonote  equilibrium  and  the  parameters  hove 
been  measured  in  the  laboratory,  all  by  means  of  the  resonotor  method  (9). 
The  laboratory  investigations  also  revealed  a  pH-dependent  relaxation 
involving  magnesium  carbonate,  which  has  been  found  to  play  a  minor  but 
significant  role  in  sea  water  absorption  ( 10). 

An  absorption  formula,  based  solely  on  known  chemical  processes  would 
be  for  too  complex  and  the  accuracy  would  be  limited  os  well.  However, 
since  the  range  of  environmental  parameters  r-  -he  World  Ocean  is  very 
limited,  simplifying  approximations  con  be  moot,  non.aly,  that  the  losses 
for  both  pH-dependent  relaxations  increase  exponentially  with  pH  and  the 
relaxation  frequencies  increase  exponentially  with  temperature.  Thorp  s 
formula  con  then  simply  be  modified  by  adding  the  third  relaxation  and 
including  the  required  pH  and  temperature  corrections. 

A  3-relaxation  formula  has  already  been  developed.  Predictions  based  on 
archival  pH  data  have  been  tested  against  all  the  available  sound-channel 
with  good  results  1 1 1,12). 


Absorption  Model 


A=A,(t1gS04)+A2(B(OH)3)+A3(MgC03) 
A„=  anf2f„/(f2+f2) 

a,  =  o.5xia'd(km)/20  f,  =  50xioT/6B 

a2=o.ixio<pM  8>  f2=o.9xioT/7° 

(pH-8)  *  .  <rtT/30 

a  _  n  m..  <  n' r  '  f  —  A  I  A 


Atlantic  4'C  pH  8.0 

A=0.Q07f  2  +0.1  f2/(1+f2)*0.t8f2/(62+y2) 

N. Pacific  4*C  pH  7.7  _  ~~ 

A=Q.007f2+0.05f2/(1+f2)+0.09  f2/(62+f2) 
Mediterranean  14#C  pH  8.3 
A=0.006f  2+G.26f  2/(1.42*f  2)+0,78  f2/( 122+f 2) 

Red  Sea  22#C  pH  8.2 

A=0-0Q4f  2*0.27f  2/(  1 .82+f  2)  ♦  1.1f2/(242+f2) 

sub- Arctic -PC  pH  8.3 

A^O.Ol  f2  +  0.17  f2/(0.852*f2)*0.24  f  2/(42*f  2  ) 

Figure  1:  Simplified  absorption  formulae. 


In  the  3-relaxation  formula  of  Figure  1,  A  is  in  dB/km,  frequency  f  and 
relaxation  frequencies  fn  are  in  kHz,  temperature  T  is  in  °C  and  pH=8.0  is 
the  reference  value.  The  pure  water  term  is  neglected,  making  the  formula 
valid  for  frequencies  less  than  roughly  100  kHz.  PH  values  in  the  World 
Ocean  vary  roughly  from  7.7  to  8.3,  which  corresponds  to  an  absorption 
rotio  of  os  much  as  4/1  at  the  lower  frequencies. 

The  magnesium  sulfate  term  includes  the  depth  factor  D(km),  which  is 
adopted  from  the  pressure  correction  of  Fisher  end  Simmons  (12).  Depth 
dependencies  of  the  other  two  relaxations  are  not  yet  known,  however, 
measurements  indicate  that  boric  acid  effects  are  negligible.  Magnesium 
carbonate  depth  effects  can  be  neglected  because  its  contribution  is  so 
small.  Salinity  dependence  will  be  considered  later.. 

Specific  coefficients  for  several  experimental  areas,  are  shown  in  the 
bottom  box.  Note  that  the  magnesium  sulfate  terms  are  approximations 
valid  only  for  frequencies  much  less  than  the  relaxation  frequencies  f,. 
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Figure  2:  Thorp's  data  and  3-relaxation  model. 

Figure  2  shows  Thorp's  data  compared  to  the  3-relaxation  model.  The 
Individual  components  are  identified  end  the  top  curve  is  their  sum.  The 
overall  fit  to  the  data  is  seen  to  be  as  good  or  better  than  with  Thorp's 
2-component  formula. 

Note  that  the  boric  acid  (A2)  coefficient  is  lower  than  that  in  the  Thorp 
formula  by  some  10$;  l.e.,  the  coefficient  in  dB/km  becomes  equal  to  the 
value  in  Thorp's  equation,  which  was  given  in  dB/kyd.  Differences  in  the 
tote 1  absorption  at  lower  frequencies  are  then  made  up  by  the  magnesium 
carbonate  component  (A3). 

The  parameter  adjustment  is  mainiy  justified  on  the  basis  of  data-fit, 
the  third  component  being  essential  to  the  model.  When  the  sea-water 
resonator  data  were  fitted  with  8  2-relaxation  model,  there  were  serious 
discrepancies  at  higher  pH  values.  Correction  of  Thorp's  equation  for  pH8.5 
gave  values  that  were  much  too  low  at  the  lower  frequencies,  which  was 
clear  evidence  of  the  existence  of  a  third  component.  Sea  water  synthesis 
experiments  were  then  carried  out.  The  mechanism  was  identified  as  a 
magnesium  carbonate  relaxation  and  the  relaxation  parameters.  Including 
temperature  and  pH  dependence,  were  measured. 
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Figure  3:  Model  and  data  comparison. 

Figure  3  compares  the  3-component  relaxation  model  predictions  and 
data  from  sound-channel  experiments  In  the  Mediterranean,  North  Pacific 
and  North  Atlantic  (Thorp). 

In  the  North  Pacific  case  1141,  the  lower  value  pHs7.7  reduces  both  the 
boric  acid  (A2)  and  the  magnesium  carbonate  (A3)  coefficients  by  a  factor 
of  two  compared  to  the  N.  Atlantic.  Relaxation  frequency  depends  only  on 
temperature  and  remains  the  same. 

In  the  Mediterranean  case  115),  the  higher  value  pHz0.3  increases  both 
the  boric  acid  (A2)  and  the  magnesium  carbonate  \A3)  coefficients  by  a 
factor  of  two  compared  to  the  N.  Atlantic.  However,  since  the  relaxation 
frequency  is  higher,  the  curves  do  not  differ  by  as  large  a  large  factor  at 
the  lower  frequencies. 

The  value  pH  8.0  tjas  been  assumed  for  Thorp's  experiment  and  is  used  as 
reference  V8lue  for  the  3-relaxatton  model.  Predictions  based  on  archival 
pH  values  then  show  agreement  within  experimental  limits  for  all  regions 
examined.  Although  no  discrepencies  are  now  apparent,  small  adjustments 
of  any  of  the  parameters  can  be  made  whenever  new  absorption  and  pH  data 
indicate  the  need. 
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PH  Profiles 


Typical  pH  profiles  for  the  North  Atlantic  and  North  Pacific  Oceans  and 
the  Mediterranean  Sea.obtained  from  the  6EQSECS  expedition  reports  1 161, 
are  shown  in  Figure  4. 

The  most  striking  features  are  the  relative  constancy  of  the  pH  values  at 
the  surface  and  the  slow  variations  at  great  depths.  Variability  with  depth 
means  that  net  absorption  realized  for  the  various  propagation  modes  will 
also  depend  very  strongly  on  the  pH  profile  and  the  particular  ray  paths 
involved. 

In  the  analysis  of  sound-channel  propagation  experiments,  axial  values 
of  pH  and  tempersture  have  proved  to  give  adequate  approximations,  within 
expected  limits  of  experimental  error.  Axial  depths  at  mid-latitudes  are 
near  I  km.  where  the  regional  variability  tends  to  be  greatest.  At  higher 
latitudes,  the  sound-channel  axis  rises  nearer  to  the  surface  where  the  pH 
values  are  uniformly  high.  For  convergence-zone  propagation,  integration 
of  loss  over  the  ray  paths  by  numerical  methods  has  shown  that  the  2  km 
depth  values  give  accurate  estimates  except  for  possible  errors  8t  higher 
latitudes  Surface  values  ere  obviously  appropriate  for  the  surface  duct. 
The  interim  global  model  has  been  based  on  these  approximations. 
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Global  Model 


A=A1(MgS04)+A2(B(0H)j)+AJ(MgC0j) 
An  =  (s/35)  a„f  2fn/(  f  2+  f 2 ) 

a,  =  o.5xia'D(‘“)/20  f,=50xioT/6° 


a2=  o.i  k 

a  3=  0.03  K 


f2  =  0.9x10 
f  3=  4.5x10 


T/70 


T/30 


Figure  5:  Global  model  absorption  formula. 


Variability  of  pH  is  clearly  the  major  limiting  factor  in  the  accuracy  of 
the  absorption  model.  In  the  proposed  global  model  of  Figure  5,  the  pH 
parameter  KslO^8*  has  been  substituted  in  the  absorption  formula  of 
Figure  I.  Salinity  dependence  is  taken  as  S/35  with  the  caveats  noted  (11); 
i.e.  errors  probably  become  excessive  outside  the  range  30-40  ppt  since 
the  changes  in  relaxation  frequency  must  then  be  considered.  This  could 
require  analysis  of  constituents,  which  is  beyond  the  present  scope. 

The  interim  model  is  based  on  estimation  of  the  effective  K  values  for 
specific  propagation  modes.  From  the  earlier  analysis,  axial  values  of  K 
and  temperature  T  (°C)  are  appropriate  for  the  sound-channel  mode  and  2 
km  depth  values  for  C2  and  other  deep  modes.  Errors  at  high  latitudes  can 
be  reduced  by  interpolation,  depending  on  the  ray  paths.  Temperature  is  not 
a  problem  since  it  can  be  derived  from  the  SVP  or  X8T  data. 

Contours  of  pH  for  the  surface  and  for  depths  0  5  km  and  1  km  8re  shown 
in  the  World  Ocean  Atlas  of  Gorshkov  ( 1 7).  Contours  for  2  km  depth  are  also 
included  in  Vol  2.  The  contour  intervals  are  0  1  pH  unit  and  interpolation 
is  required.  Correction  to  in-siitt  pressure  is  also  necessary. 

The  GEOSECS  d8la  have  also  been  used  to  estimate  the  2  km  contours  in 
parts  of  the  Pacific  Ocean  not  covered  by  the  Russian  work.  Oiscrepencies 
between  the  data  sets  have  been  arbitrarily  resolved  by  adjustment  of  the 
contours  to  minimize  effects  on  estimation  error. 

The  sound-channel  K  contour  chart  of  Figure  6  is  based  on  the  analysis 
of  Russian  pH  contours  by  Lovett  (18).  Th8  CZ  and  surface  charts,  derived 
mainly  from  the  Russian  report,  8re  shown  in  Figures  7  and  8,  respectively 
The  CZ  chart  is  a  modification  of  the  earlier  2  km  chart  with  corrections 
for  high-latitude  effects.  Effective  K  values  were  obtained  by  integration 
over  appropriate  ray  paths 


-  6  - 


The  absorption  formula  of  Figure  5  has  been  based  on  the  results  of  both 
laboratory  and  sound-channel  experiments.  Because  the  frequency  range  of 
the  resonator  was  limited  to  10-100  kHz,  it  was  not  possible  to  measure 
the  boric  acid  relaxation  parameters  by  this  means.  The  field  dote  were 
primarily  used  in  determining  the  absorption  spectra  at  lower  frequencies. 
Thorp's  date  of  figure  2,  considered  to  be  the  most  accurate,  was  heavily 
weighted  in  the  analysis. 

The  pH  data  for  the  sound-channel  experiments  was  archival  and  errors 
are  estimated  to  be  of  the  order  ±0.05  units,  which  corresponds  to  ±12* 
error  in  magnitude  at  lower  frequencies.  Therefore,  the  absolute  accuracy 
of  the  pH-dependeni  pert  of  the  formula  is  expected  to  be  of  this  order. 

The  accuracy  of  the  pH-independent  part  is  as  good  or  better.  Resonator 
studies  of  sea-water  absorption  in  the  low  kHz  range  have  been  carried  out 
recently  at  the  Academia  Sinica  laboratories  and  measurements  over  the 
pH  range  of  interest  confirm  overall  formula  predictions  with  equivalent 
accuracy  (19,201. 

Accuracy  of  field  predictions  will,  in  general,  be  limited  by  the  pH  data. 
The  K  contours  of  Figures  6-8  should  provide  reasonably  good  estimates  of 
absolute  pH  effects  for  the  three  propagation  modes  throughout  most  of 
the  World  Ocean.  Possible  exceptions  ore  regions  where  some  significant 
discrepancies  between  the  Russian  and  6E0SEC5  values  have  been  noted 
end  also  regions  where  pH  profile  changes  are  rapid.  Questions  can  also 
arise  when  there  is  no  clearly  dominant  propagation  mode. 

While  the  global  model  can  certainly  account  for  gross  regional  changes 
in  absorption,  the  overall  accuracy  may  be  inadequate  in  some  coses.  For 
example,  it  can  be  difficult  to  determine  what  values  to  use  for  arbitrary 
ray  paths.  Interpolation  between  the  sound-channel,  convergence-zone  and 
surface  values  should  help  to  improve  matters;  however,  the  errors  could 
easily  become  unacceptable  when  there  are  rapid  changes  over  the  depth 
range  in  question.  Some  subjectivity  and  uncertainty  is  always  involved  in 
such  a  process.  Furthermore,  the  errors  ere  additive  and  combined  effects 
of  ray-path  approximations  ond  errors  in  values  could  become  excessive. 
Therefore,  on  obvious  method  of  simplifying  the  computational  process,  as 
well  os  improving  the  overall  accuracy,  is  to  integrate  loss  over  oil  the 
ray  paths  using  an  appropriate  profile  for  the  region. 

tn  the  following  sections,  all  of  the  6E0SECS  dote  are  onalyzed  for  the 
purpose  of  assessing  the  occuracy  of  the  K  contour  charts.  The  second  and 
perhaps  major  purpose  is  to  investigate  K-profiie  estimation  methods  for 
incorporation  in  existing  computer  programs  for  propagation  loss. 
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GEOSECS  Data  Analusis 


The  Geochemical  Ocean  Sections  Study  (GEOSECS)  [161  atlas  volumes 
contain  records  of  the  oceanographic  data  token  during  the  International 
Decade  of  Ocean  Exploration  (1DOE)  1970-1980.  The  areas  covered  include 
the  Atlantic,  Pacific  and  Indian  Oceans.  The  expedition  tracks  followed 
depth  contours  of  4  km  or  more  and  coverage  is  not  comprehensive.  Data 
are  missing  for  many  of  the  stations  also,  particularly  the  first  part  of 
the  North  Atlantic  track. 

The  Carbonate  Chemistry  sections  of  each  volume  contain  the  tobies  of 
temperature,  salinity,  pH  and  associated  parameters  vs  depth.  The  pH  data 
is  the  primary  concern  in  this  work.  Procedures  and  method  of  pH  analysis 
are  discussed  in  the  atlases  and  references.  Values  corrected  to  in-situ 
temperature  and  pressure  ore  included  and  ore  used  in  this  report. 

Data  for  the  three  oceans  are  presented  in  sections,  each  beginning  with 
a  track  chart  showing  the  station  numbers.  Station  sequence  is  followed 
except  for  the  Indian  Ocean  section  where  the  profiles  404-407,  tor  the 
Mediterranean  Sea,  the  Red  Sea  and  the  6ulf  of  Aden,  are  placed  at  the  end 
in  reversed  order. 

The  graphs  at  the  *op  of  the  figures  following  contain  5  sequential  pH  vs 
depth  Profiles.  Depth  intervals  were  chosen  to  give  the  minimum  required 
detail  and  the  original  daio  were  interpolated  to  obtain  the  values  when 
necessary.  Any  fine  structure  present  was  minimized  by  smoothing.  The  pH 
scale  for  each  station  is  indicated  by  solid  vertical  lines  with  the  value 
pHG.O  at  the  top.  The  tick  marks  are  0.1  pH  units  and  successive  profiles 
are  displaced  two  units  to  the  right.  The  tables  below  the  graphs  show  the 
station  numbers,  station  coordinates  end  calculated  K  values  vs  depth. 

The  K  contours  at  five  selected  oept'  s  0,  0.5,  I,  2  end  4  km  ore  shown  in 
Figures  41-45  Tho  contours  were  derived  by  plotting  the  K  values  on  the 
charts  and  adiusting  for  'best-fit*,  taking  into  consideration  values  and 
trends  indicated  by  the  Russian  profits.  Note  the*,  the  K  values  refer  to 
regions  separated  by  contour  lines  md  not  the  contour  lines  themselves, 
os  in  the  earlier  cherts.  This  simplifies  selection  of  values  in  the  regions 
of  gradual  change  and  serves  to  identity  problem  regions  of  rapid  change 
as  well. 

The  five  K  values  provided  by  the  charts  should  permit  quite  accurate 
estimation  of  the  K  profiles  in  regions  whore  no  data  ore  available.  In 
later  sections,  several  alternative  methods  for  formulating  profiles  ere 
proposed.  The  expected  error  in  predicted  loss  due  to  profile  estimation 
errors  is  also  considered. 
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Atlantic  Ocean 


TRACK  OF  R/V  KIMQRR, 
GEOSECS  ATLANTIC  EXPEDITION,  1  972-73 


Figure  9:  Atlantic  track. 
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Figure  10:  Atlantic  Ocean  stations. 
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Figure  1 1:  Atlantic  Ocean  stations. 
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Figure  i2:  Atlantic  Ocean  stations. 
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Figure  13.  Atlantic  Ocean  stations. 
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Figure  14:  Atlantic  Ocean  stations. 
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Figure  16:  Atlantic  Ocean  stations. 
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Figure  20:  Atlantic  Ocean  stations. 
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Pacific  Ocean 


TRACK  OF  R/V  MELVILLE, 
GEOSECS  PACIFIC  EXPEDITION,  1373-7A 


Figure  21:  Pacific  track. 
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Figure  24:  Pacific  Ocean  stations. 
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Figure  25:  Pacific  Ocean  stations. 
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Figure  26:  Pacific  Ocean  stations. 
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Figure  27:  Pacific  Ocean  stations. 
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Figure  28:  Pacific  Ocean  stations. 
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Figure  29;  Pacific  Ocean  stations 
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Figure  30;  Pacific  Ocean  stations. 
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Figure  31:  Pacific  Ocean  stations. 
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Figure  32:  Pacific  Ocean  stations. 
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Figure  33:  Indian  Ocean  trock. 
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Figure  34:  Indian  Oc 


Figure  35;  Indian  Ocean  stations. 
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Figure  37:  Indian  Ocean  stations. 
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Figure  39:  Indian  Ocean  stations. 


Figure  43:  GEOSECS  I  km  K  contours. 
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K  Profiles 


Comparison  of  the  K  contours  indicate  no  serious  discrepancies  between 
the  new  and  earlier  versions.  However,  adjustments  of  values  end  contours 
to  fit  the  GEOSECS  data  should  improve  accuracy.  Considerable  uncertainty 
was  involved  in  developing  the  earlier  2  km  contours  because  the  Pacific 
latitudes  165°E-130',W  were  not  covered  by  Russian  charts.  Interpolation 
and  pressure  corrections  were  involved  also.  Therefore  GEOSECS  surface 
and  2  km  charts  should  be  the  more  accurate. 

The  K  contours  should  also  allow  reasonably  accurate  of  the  profiles  in 
most  of  the  World  Ocean.  A  graphic  method  \o  to  connect  the  five  points 
with  a  fair  curve,  which  can  be  quite  subjective.  An  automated  method  has 
been  developed  using  on  algorithm  suggested  by  Dr.  A.  H.  Nuttall,  which  has 
proved  to  give  good  results. 

The  five  K  points  for  depths  D=0,  0.5,  1,2,  and  4  krn  are  numbered  0-4 
and  the  profile  is  fitted  by  the  expression: 

K(D)=K(4)+  (C0  ♦  Ct  D  ♦  C2  D2  ♦  Cj  D5  ♦  C4  D4]  expHa  D)b] 

where  8=1 /km  and  b=1.5  have  Lsen  chosen  by  trial  to  obtain  ‘best"  results. 

The  five  equations  for  n=0,l, 2,3,4  to  be  solved  are  then  given  by 

CQ  ♦  C,  Dn  ♦  C2  Dft  *  ♦  C3  Dft 5  ♦  C4  D/  *  lK(Dn)  -  K(4))  exptDj  3) 

and  solution  for  the  coefficients  con  be  obtained  either  algebraically  or  by 
writing  the  equations  in  matrix  form  end  inverting. 

Figures  46-49  show  some  results  for  the  North  Atlantic,  North  Pacific, 
Indian  and  Southern  Oceans.  The  circles  indicate  the  curve-fit  to  the  five 
chart  values  for  each  region.  A  computer  search  is  then  made  of  the  entire 
library  to  find  the  five  GEOSECS  profiles  with  least-square  errors  for  the 
given  points.  The  error  6  rms  is  shown  for  each  cose  Selected  profiles  ore 
identified  by  station  number  onu  error  Since  the  error  only  indicates  the 
goodness  of  fit  at  the  five  points,  actual  errors  for  particular  ray  paths 
con  obviously  be  greater  or  less 

For  the  all  the  cases  shown,  the  selected  profiles  fall  within  the  regions 
in  question  and  the  errors  appear  to  be  minimal  In  other  cases  where  the 
method  recovers  profiles  outside  the  area,  the  errors  ore  equally  smoll  An 
alternative  procedure  would  be  to  select  the  profile  with  the  best  overall 
fit  in  place  of  the  algorithm  profile. 
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Figure  48:  Indian  Ocean  K  profiles. 


Figure  49:  Southern  Ocean  K  profiles. 
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Temperature  Profiles 


Figure  50:  Typical  temperature  profiles. 

Temperature  profiles  are  also  required  for  the  ray-integration  method. 
Figure  50  shows  some  typical  examples  for  the  regions  to  be  considered. 
Temperature  profiles  are  obviously  quite  similar  to  K  profiles  and  can  be 
formulated  by  the  same  algorithm  method.  Salinity  variations  with  depth 
are  small  and  can  be  Ignored. 

Figures  51-53  show  comparison  of  the  temperature  and  K  profiles  for  all 
reported  stations  in  the  three  oceans.  Station  sequences  read  from  left  to 
right  and  are  the  S8me  as  before.  Starting  and  ending  station  numbers  are 
indicated  on  each  graph.  Note  that  the  spacing  has  been  increased  In  the 
lower  Indian  Ocean  graph  to  eliminate  overlap  of  the  Gulf  of  Aden,  Red  Sea 
and  Mediterranean  Sea  temperature  profiles,  which  might  otherwise  cause 
confusion. 

While  there  is  a  considerable  degree  of  similarity  between  the  profiles 
In  many  regions,  the  overall  correlation  is  generally  poor  and  clearly  not 
sufficient  to  permit  estimation  of  K  profiles  from  temperature  profiles 
with  the  required  accuracy. 
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Figure  51:  Temperature  and  K  profile  comparison. 
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Figure  53:  Temperature  and  K  profile  comparison. 
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Model  Comparison 


Absorption  losses  (dB)  calculated  by  both  the  K  model  and  Thorp  equation 
for  sound-channel,  convergence  and  surface-duct  propagation  modes  are 
shown  in  Figures  54-56.  The  differences  A=A(Kmod)-A(Thorp)  include  the 
expected  errors  of  the  K  model,  indicating  the  degree  of  significance. 

The  K  model  calculations  used  ray  integration  and  the  algorithm  profile. 
Values  of  K  and  temperature  at  the  five  depths  were  fitted  to  generate  the 
profiles.  Salinity  was  taken  os  constant  S=35  except  for  the  Mediterranean 
cose  where  S=38  was  used.  Absorption  loss  at  selected  frequencies  F(kHz) 
at  ranges  R(km)  was  calculated  by  integration  over  appropriate  ray  paths. 
Expected  errors  were  calculated  concurrently  using  AK=+0.05  for  the  sound 
channels  and  CZ  and  AK=i0.1  for  the  surface  ducts.  Thorp  values  were  also 
calculated  concurrently  to  minimize  relative  errors. 

Figure  54  shows  one-way  sound-channel  losses  at  R=500  km.  The  large 
range  of  A  values  reflects  the  regional  pH  variability  at  axial  depths. 

Figure  55  shows  two-way  absorption  losses  for  active  sonar  in  the  CZ 
mode.  Since  the  major  part  of  the  ray  paths  falls  well  below  the  depths  of 
high  pH  variability,  A  values  tend  to  be  small  in  both  magnitude  and  range. 

Figure  56  shows  two-way  surface-duct  losses  at  R=50  km.  The  A  values 
are  uniformly  high  because  of  the  high  surface  pH. 
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Figure  54:  One-way  sound-channel  absorption  values. 
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Figure  55:  Two-way  convergence-zone  absorption  values. 
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Figure  56:  Two-way  surf8ce-duct  absorption  values. 


Conclusions  and  Recommendations 


Two  methods  for  estimating  pH  effects  on  sound  absorption  have  been 
proposed  for  use  in  conjunction  with  the  absorption  formula  of  Figure  5. 
The  interim  method  employs  K  contours  for  specific  propagation  modes; 
i.e.  sound  channel,  CZ  and  surface  duct  (Figures  6-8).  This  method  is  useful 
for  rapid  calculation  but  involves  judgement  if  there  is  no  single  dominant 
mode.  The  second  method  employs  K  profiles  constructed  from  K-contour 
charts  at  five  depths  0,  0.5,  t,  2  and  4  km  (Figures  41-45).  This  method  is 
useful  for  integrating  loss  over  ell  ray  paths  by  computer  methods. 

One  purpose  of  the  GEOSECS  data  analysis  has  been  to  investigate  the 
consistency  of  the  pH  data.  Results  indicate  agreement  between  the  new  K 
contours  and  the  earlier  versions.  The  sound-channel  contours  ore  in  good 
agreement  with  Lovett's  charts.  Agreement  with  Russian  charts  at  fixed 
depths  is  also  good  except  for  possible  minor  discrepencies  in  surface 
values,  which  may  reflect  seasonal  variability.  6E0SECS  data  were  taken 
during  warmer  periods  and  surface  temperatures  appear  to  be  consistently 
high,  while  the  Russian  data  are  seasonally  averaged.  Within  these  limits, 
errors  for  both  surface  charts  appear  comparable. 

The  second  purpose  has  been  to  devise  a  simple  procedure  for  estimating 
the  K  profiles  required  for  propagation  analysis  by  numerical  methods.  By 
integration  losses  over  all  ray  paths  using  a  profile,  all  uncertainty  about 
what  values  to  use  is  eliminated.  The  five  points  from  the  contour  charts 
should  allow  reasonably  accurate  approximation  of  K  profiles  throughout 
most  of  the  World  Ocean  in  regions  where  no  ectual  data  can  be  found. 

Profiles  con  be  generated  graphically  by  drawing  a  fair  curve  through  the 
five  points.  Since  making  the  fit  consistent  with  actual  profiles  involves 
some  subjectivity,  the  algorithm  method  may  be  preferred. 

An  alternative  method  is  to  do  a  computer  search  a  library  of  profiles  and 
select  the  one  with  least-square  error  for  the  five  points  in  question.  This 
procedure  will  usually  recover  a  profile  in  the  some  generol  area.  If  not, 
discrepencies  between  known  and  recovered  profiles  ere  usually  found  to 
be  minimal.  Another  alternative  is  to  check  the  overall  fit  of  the  recovered 
profiles  and  arbitrarily  pick  the  one  that  appears  best  suited  to  the  region 
under  consideration. 

One  of  the  main  advantages  of  the  library  method  is  that  additions  can  be 
mode  as  new  data  become  available  While  some  revisions  of  the  contours 
may  then  be  required,  the  extrapolation  errors  involved  in  approximating 
profiles  should  be  reduced.  This  method  also  seems  to  be  suitable  for  easy 
implementation  in  existing  computer  programs. 
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